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INTRODUCTION 

This  paper  is  about  finite  element  calculations  of  the  Mode  1 
stress  intensity  factor  (ifj),  and  the  T-stress  (J),  for  single 

edge-notched  tension  (SENT)  specimens.  Linear  elastic 
fracture  mechanics  (LEFM)  is  used,  peases  of  traction 
application^- method  and  plane  state  (TAMPS)  are 
considered,  applied  uniform  load  and  plane  stress  (LS), 
applied  uniform  load  and  plane  strain  (LN),  applied  uniform 
displacement  and  plane  stress  (DS),  applied  uniform 
displacement  and  plane  strain  (DN).  For  all  cases,  the 
material  is  considered  to  be  incompressible,  i.e.,  0.5  is  used 
for  Poisson's  ratio  (v).  Young's  modulus  is  referred  to  as  E. 
The  total  height  (perpendicular  to  the  crack)  of  a  SENT 
specimen  is  called  h ,  and  the  total  width  (parallel  to  the 
crack)  is  called  w .  The  crack  length  is  called  a. 


FINITE  ELEMENT  CALCULATIONS 

The  ABAQUS  finite  element  program,  version  5.6  was  used. 
All  calculations  were  for  a  linear  material  with  small  strains 
and  small  displacements.  Figure  0  shows  the  general  form  of 
the  finite  element  meshes.  One  half  of  each  SENT  specimen 
was  modeled,  cut  through  the  line  of  symmetry,  the  crack 
line.  In  Figure  0,  the  'y's  along  part  of  the  bottom  of  the 
mesh,  indicate  that  each  point  ahead  of  the  crack  tip  along  a 
mesh's  bottom  edge  was  constrained  from  displacement  in 
the  y  direction.  The  'x's  along  the  top  of  the  mesh  indicate 
that  each  point  along  a  mesh's  top  edge  was  constrained 
from  displacement  in  the  x  direction.  (The  top  edge  x 
constraint  was  used  to  simulate  the  gluing  of  an 
experimental  specimen  to  an  end  fixture.)  Each  mesh  was 
divided  into  16  element  sectors  of  approximately  equal  angle 
radiating  from  the  crack  tip.  Each  mesh  was  divided  into  64 
element  orbits  about  the  crack  tip.  Each  mesh  consisted  of 
16*64  =  1024  elements.  The  elements  were  rectangular  8 
node  (one  node  at  each  corner,  one  node  on  each  side) 
isoparametric.  Each  crack  tip  element  was  formed  by 
positioning  the(S/iodes  along  one  side  at  the  crack  tip.  Ail  of 
the  crack  tip  nodes  were  tied  together.  Element  corner 
nodes  were  located  at  the  intersection  of  the  sector  lines  and 
orbit  curves.  Element  side  nodes  were  located  midway 
between  the  corner  nodes.  The  nodes  on  the  radial  sides  of 
the  crack  tip  elements  were  located  at  the  quarter  points. 


The  family  of  o.oif^^in^d^nconstructed  so  that  the 
innermost  curve  was  circular.  The  shapes  of  the  curves 
gradually  became  more  rectangular  nearer  to  the  mesh 
boundaries  (the  two  sides  and  the  top).  Beginning  at  the 
crack  tip,  the  distances  between  adjacent  orbit  curves 
increased  linearly  in  the  -x ,  -kr,  and  +y  directions.  The 

radius  of  the  innermost  (circular)  orbit  curve  was  wj 104  .(396^ 
(4  TAMPS  *  11  h/w  *  9  a/w)  ABAQUS  problems  were 
solved  to  generate  the  included  data.  A  FORTRAN  90 
program  was  written  and  used  to  automatically  generate  a 
mesh  and  ABAQUS  input  file  from  the  values  of  a  number  of 
variables  in  a  parameter  file.  ABAQUS  was  run  on  a  UNIX 
machine.  Typically,  99  runs  were  made  at  one  time  from  a 
batch  file  running  overnight. 

THEORY  AND  NUMERICAL  IMPLEMENTATION 
Figures  1-4  show  the  T-stress  parameter/}  plotted  versus 
the  crack  ratio.  /}  is  defined  as 

Tjm 

p  kT 

Kx  is  equal  to 

J  is  the  energy  release  rate  at  the  crack  tip.  E'  is  defined 
as 

E  (plane  stress) 

— (plane  strain)  *  ^ 
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ABAQUS  has  a  command  to  calculate  J,  and  uses  the 
energy  domain  integral  method.  For  each  ABAQUS  run,  64 
calculations  of  J  were  made,  one  for  each  element  orbit 
curve.  For  each  run,  all  64  J  values  were  close  to  identical. 

Williams  [1]  showed  that  the  stresses  in  a  linear  elastic  body 
can  be  written  in  an  infinite  series  expansion.  Near  the  crack 
tip  the  effect  of  the  higher  order  terms  can  be  neglected,  and 
can  be  written  as 
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Along  the  crack  flank  ( 6=  n),  the  first  term  in  Equation  (4)  is 
zero,  and  becomes 

q^  =  r.  (5) 

Therefore,  along  the  crack  flank,  and  near  the  crack  tip,  a 
finite  element  calculation  of  oa  should  be  approximately 
equal  to  T.  Ayatollah!  [2]  has  shown  that  a  better 
approximation  for  T  can  be  determined  by  using 

a  =T  =  E'^,  (6) 

dx 

where  ux  is  the  displacement  along  the  crack  flank  in  the  x 
direction.  For  the  data  of  this  paper,  T  was  calculated  using 
Equation  (6).  dux/dx  was  approximated  by  fitting  a  least 

squares  line  to  the  data  points  (x,ax)  of  the  nodes  on  the 
crack  flank  within  the  distance  of  a/10  of  the  crack  tip. 

Figures  5-8  show  normalized  T-stress  [T/a^  )  plotted  versus 
the  crack  ratio.  is  the  specimens  average  end  stress; 

i.e.,  the  sum  of  the  mesh  top  edge  nodal  forces  in  the  y 

direction  divided  by  w .  (For  cases  LS  and  LN  these  nodal 
forces  were  applied,  and  for  cases  DS  and  DN  they  were 
induced.) 

NUMERICAL  RESULTS 

Each  of  Figures  1-8  shows  11  curves,  one  for  each  of  the 
specimen  aspect  ratios  h/w-{  1/4,  1/3,  1/2,  1/1,  2/1,  3/1, 
4/1,  5/1,  6/1,  7/1,  8/1}.  Each  curve  was  constructed  using  9 
points,  for  the  crack  ratios  a/w -{0,05,  0.10,  0.20,  0.30, 
0.40,  0.50,  0.60,  0.70,  0.80}. 

Notice  that  for  a  SENT  experimental  specimen,  both  T  and 
Kx  can  be  determined  using  the  figures.  Since  ol  is  either 
applied  or  measured,  one  of  Figures  5-8  will  give  T.  Once 
T  is  known,  the  corresponding  figure  of  Figures  1-4  will  give 
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Figures  1  and  5  were  compared  to  Figure  4  from  the  T- 
stress  compendium  by  Sherry  [3],  and  showed  reasonable 
agreement. 

The  figures  generally  show  that  the  fundamental  characters 
of  the  curves  change  when  the  specimen  aspect  ratio  (h/w) 

is  less  than  1.  There  is  little  difference  between  the  effects  of 
plane  stress  and  plane  strain  when  uniform  load  is 
prescribed  (cases  LS  and  LN).  However  when  uniform 
displacement  is  prescribed  (cases  DS  and  DN)  the 
differences  between  plane  stress  and  plane  strain  can  be 
significant,  especially  for  h/w<  1.  The  figures  show  that 
both  f}  and  T/aM  crucially  depend  on  whether  uniform  load 
or  uniform  displacement  is  applied  to  a  specimen.  One 
condition  should  never  be  used  to  approximate  the  other. 


CONCLUSIONS 

The  boundary  condition  of  Figure  0  (in  which  the  entire  mesh 
top  edge  was  constrained  in  the  x  direction)  was  used 
because  it  seems  to  be  the  only  condition  that  can  be 
practically  implemented.  Devising  an  experimental  method  to 
constrain  only  one  point  in  the  x  direction  seems  difficult. 
The  result  is  that  for  values  of  h/w  <  1,  relatively  huge  T- 
stresses  can  be  generated. 

Applying  a  load  to  a  SENT  specimen  that  is  even 
approximately  uniform  seems  to  be  an  impossibility,  unless 
the  specimen  is  extremely  long.  Even  if  the  end  fixtures  are 
pinned  at  their  middles,  the  most  that  can  be  said  is  that  the 
centroid  of  the  end  loads  are  located  at  the  pin  centers. 

For  experimental  single  edge-notched  tension  specimens, 
calculations  using  uniform  applied  load  should  never  be 
substituted  for  calculations  using  uniform  applied 
displacement  when  determinining  values  of  the  stress 
intensity  factor  (Kx)  and  the  T-stress  (r). 
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Fig.  0  Typical  Finite  Element  Mesh. 
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